Asymmetric fluctuation- relaxation paths in FPU models 
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Abstract 

\ A recent theory by Bertini, De Sole, Gabrielli, Jona-Lasinio and Landim predicts a temporal 

asymmetry in the fluctuation-relaxation paths of certain observables of nonequilibrium systems 
' in local thermodynamic equilibrium. We find temporal asymmetries in the fluctuation-relaxation 

o , 

, CJ, . paths of a form of local heat flow, in the nonequilibrium FPU-/? model of Lepri, Livi and Politi. 
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I. INTRODUCTION 

In Refs. , Bertini, De Sole, Gabrielli, Jona-Lasinio and Landim proposed a gener- 
alization of Onsager-Machlup's theory 3] to the large fluctuations around nonequilibrium 
steady states. This theory predicts temporal asymmetries between fluctuation and relax- 
ation paths of certain observables, and includes the results of Derrida, Lebowitz and Speer 
a. a.d of Bodineau and De„ida .I as speca. case. In pa.*u.a... the theory of J Q 
concerns stochastic lattice gases, which admit the hydrodynamic description 



T^{q) , Q= Q{u,t) , (1) 



where g is the vector of macroscopic observables, u and t are the macroscopic space and 
time variables, D is the Onsager diffusion matrix, x is the linear response conductivity, and 
V-E is the external driving force. For such models, Refs. prove that the spontaneous 
fluctuations out of a steady state most likely follow an adjoint hydrodynamic equation: 

dtg = V*{g) , with V*{g)=V{g)-2A, (2) 

where T) has been decomposed as 

'D{g) = \v- (x{Q)'^^-f] (3) 



2 V J 

and ^ is a vector field orthogonal to the thermodynamic force, i.e. to the functional derivative 
of the entropy with respect to the state 6S/6g. The fact that A is orthogonal to the 
thermodynamic force implies that it does not contribute to the entropy production: A is 



the non-dissipative part of _t 
Because the theory of 



le dynamics. 



is developed for stochastic processes, the question arises 
whether its predictions can be observed in the dynamics of time reversal invariant, deter- 
ministic, particle systems, such as those of nonequilibrium molecular dynamics (NEMD) 
0]. Indeed, the stochastic description of a given system is often thought to be a reduced 
representation of the microscopic deterministic dynamics of its particles. One difficulty in 
addressing this question is that not all reductions of a microscopic dynamics to a meso- 
scopic (stochastic) one are equivalent, and it is not always easy to identify the stochastic 
counterparts of the observables of deterministic systems. 

A nonequilbrium steady state requires at least two reservoirs at different thermodynamic 
states, in contact with the system of interest, and NEMD provides techniques to avoid the 



simulations of the reservoirs. However, present day NEMD simulations cannot be performed 
with sufficiently many particles, and for sufficiently long times to achieve the mesoscopic 
level of description. Therefore, in Ref. [3| the simplest of particle systems, the nonequilibrium 
Lorentz gas, was investigated, but no temporal asymmetry was found in the fluctuation- 
relaxation paths (FRPs) of its current. Here, we consider the FRPs of the local heat flux 

nn 

of a more realistic model: the nonequilibrium FPU-/? model of Lepri, Livi and Politi |8|, |9||, 
and we observe that they are asymmetric in time. The model consists of the usual FPU-/? 
chain with anharmonic oscillators of equal masses m, located at Xj, j = 1,...,N, and 
with two Nose-Hoover "thermostats" at different temperatures acting on the first and the 
last oscillator. The interaction potential and the internal energy are thus given by 

2 4 ^1 

Viq) = \ + P\; and H = ^^nq] + V{q,^, - q,) , (4) 

j=0 

where qj = xj — ja is the displacement of the j-th oscillator from its equilibrium position, 
and a is the equilibrium distance between two nearest neighbours. The equations of motion 
take the form 

qi = F{qi - go) - F{q2 - qi) - OAi , qN = F{qN - qN~i) - F{qN+i - qN) - (rqN (5) 
(ij = F{qj - qj^i) - F{qj+i - qj) , for j = 2, - 1 (6) 

where F{q) = —V'{q), and the terms (iXr obey 



with the boundary conditions go = qN+i = 0. We take (3 = 0.1 and m = 6r = 0^ = 1, 
as in .8]- Equations (jSj) and ((Zj) constitute the Nose-Hoover thermostats at "temperatures" 
Tr and T^, with response times^ and 6i. This nonequilbrium FPU-/5 model is dissipative 
[ll[,but time reversal invariant |l2l|. Note that in our context (T^ — Ti) is just a parameter 
expressing the distance from the equilibrium state, and not a real temperature gradient. 

Let L be a number of contiguous oscillators in a sublattice L = {ji, . . . , jl} of the chain, 
whose center coincides with the center of the chain, and lei ii = L / N he the fraction of the 

Q 

chain occupied by L. The instantaneous local heat flux is defined by 

1 1 

= - ^ - {xj+i - Xj) (gj+i + qj)F{q,j+i - g^) + q^hj, (8) 



where 



1 1 



and -F(q') = —V'{q). There is another definition of local heat fiux, which could be considered: 

1 -^^^ 1 1 
j=ji 

The first term of / is the average over L of the heat fiux in the limit of small (compared to 
the lattice spacing a) oscillations as in ^S], while the second term gives the fiux through the 
boundary of L. The quantities and / are supposed to be equivalent in the limit of small 
oscillations and long chains. Therefore, for the sake of brevity, in this paper we do not study 
the properties of the FRPs of /, which will be considered in a future paper [l^, along with 
the FRPs of other observables. 

Because of the large difference between the mesoscopic level of description, and the micro- 
scopic level simulated by us, it is not obvious that our FRPs are directly related to those of 



1|, 12]; if they are not, their temporal asymmetry is a feature of their dynamics which, to the 
best of our knowledge, has not been studied before in microscopic deterministic dynamics. 
However, the fact that our results do not change qualitatively with the growth of N suggests 
a relation between our results and the asymmetries of the large deviations of Refs. hj, ^1 , and 

riQ 

indicates how the theory of |2j may extend to the microscopic level of description. 



II. DETERMINISTIC FLUCTUATION-RELAXATION PATHS 

Given the observable J-' : Ai —>■ IR, the identification of a FRP around its steady state 
(expected) value IE{J^) requires some care. Indeed, the time series of JF looks very noisy, 
and takes the value iE(jF) only at discrete instants of time. Similarly, the fiuctuation values 
selected to be observed are also achieved only at discrete instants of time. Averaging over 
many particles, and over many microscopic times, stabilizes the signal, but makes negligible 
the frequency of the FRPs, to the point that it becomes impossible to assess their statistical 
properties. 

By analogy with the notation of we denote by J^t the quantity J^{S^T), if F is the 
initial point of the trajectory in Ai, and : A4 ^ Ai is the time evolution operator. We 
want to identify the most likely FRP, according to the SNS probability distribution, starting 
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at Ti — TE{!F\ reaching the fluctuation value — a-nd returning to Tf — IE{T) 

after a certain time. To asses the dependence of the phenomenon on N, we consider two 
kinds of fluctuation values: 

T;{T) = 1E{T) + sa{T), and T}^\J^) = 1E{T) + -a{J^)VN; s = 1, 2, 3... (10) 

s 

with a{J^) the standard deviation of T (which depends on N), and we introduce the following 
definition of FRP: 

Definition 1. Let !F take the value T{!F) at time i, with [dj-'t/dt]^ > 0. Lett be the smallest 
time after i such that J-'t = T{!F). Choose a time tq > 0, representing the expected duration 
of the fluctuation path. The union {^i-r^ ^ [O^'^'o]} U {^t+r, t G [0)''o]} of the trajectory 
segments of duration Tq preceding t, and following t, is called an FRP. 

Numerically, the time i can be identified testing the condition 

{T{t+h) — '^){^t — ^) < 0, where h — time step, T = fiuctuation value. (11) 

Whenever a FRP is obtained, we shift it in order to have i — Given n FRPs, we 
collect them in the set {(r, J-t^'^)^ s = 1, 2 . . . , n}, and build a histogram by partitioning the 
rectangle [— To,ro] x [min^ JFJ:*\ max^ JF^''-'] with small rectangular bins, and evaluating 
the frequency of visitation of each bin. The number of temporal intervals in which [— Tq, Tq] 
is subdivided is denoted by hr, while the number of bins of \mm.T-^s ^t'\ niax,-_s .Ft*'*] is 
denoted by hj^. The area of each rectangular bin A^^^, with p = 1, . . . , 6,-, q — 1, . . . , bj^, 
is A = 2To[max,-_s J^r''^ — min^^^ V(^i-^'.f) , and the corresponding histogram is a two- 
dimensional surface, whose "crest" , defined by 

Cip) = ((2to + l)nA)-^ max « {(r, ^(^)) e A^,,, s = 1, 2 . . . , n} , p = 1, . . . , 6„ (12) 

q=l,...,bjr 

with the number of FRPs in {.}, represents the most likely FRP in the sample. To assess 
the asymmetries of the FRPs, and of the corresponding crests, we introduce the following 
definition: 

Definition 2. Let T be the fluctuation value chosen for J^. If ^{i,i,To) = {J-i-r^ 'T ^ 
[0, To]} U {J-t+T-i T G [0)'''o]} is an FRP; its asymmetry coefficient is given by 



Let C : {1, . . . br} IR he the crest obtained form the given sample of FRPs; the asymmetry 
coefficient ofC is given by 



The asymmetry coefficient a of a single FRP is defined in the same way, replacing the 
values C{p) with the corresponding values of the FRP. A FRP or a crest is symmetric if its 
asymmetry coefficient vanishes. 

Observe that the asymmetry coefficients can take positive as well as negative values, because 
the differences which define them are not taken in absolute value. This is motivated by the 
fact that random oscillations of a fluctuation path around its symmetric relaxation path 
should result in a vanishing asymmetry coefficient. 

III. NUMERICAL RESULTS 

We deal with an observable which is odd under the time reversal operation. Unfortu- 
nately, the crests obtained from our simulations are strongly noisy, as evidenced by Fig.l 
in one of the cases considered, and their asymmetry coefficient Oc is not particularly mean- 
ingful. Therefore, we study the statistics of the asymmetry coefficient a of the single FRP, 
considered as a random variable which takes different values on the different paths. Our 
results are summarized in Tables I and II, in which the fraction of oscillators over which the 
averages are computed to the total number of oscillators are fixed, but the other parameters 
defining the model vary. We conclude that the most likely FRP, which is the object of 
the theory of , is hard to assess directly, because the crest obtained for each different 
case is highly noisy and not sharp: indeed, the supports and the variances of the marginal 
probability distributions, for fixed r in the flux space, are quite wide (cf. Tables I and II). 
On the other hand, all the data Tables I and II indicate that the FR paths of have a 
positive asymmetry, like those of / A positive asymmetry of the FRPs is testified not 
only by the fact that the mean asymmetry is a positive number, but especially by the fact 
that the probability of being postive is well above 1/2. Furthermore, the dependence on the 
various parameters defining the model is rather weak. 

Therefore, we have found evidence of a new phenomenon in the nonequilibrium FPU-/5 
model of Lepri, Livi and Politi, which shows how irreversible phenomena can be obtained 




(14) 
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FIG. 1: Crest of the observable with iV = 100, r/ = 0.16, = 200, T,. = 20, ~ 24.741, and 

7^* ~ 1738. The asymmetry coefficient of the crest, computed according to (jl^f) . is ac(C) = 0.0142. 





JE{T) 


a{T) 


P(a > 0) ; Tg* 


IE{a);T* 


P{a > 0) ; ri^) 


ma) ; ^ 


25 


49.2262 


1449.5949 


0.7001 


0.0254 






50 


37.3129 


777.1178 


0.6316 


0.0336 






100 


26.4482 


569.0022 


0.6354 


0.0391 


0.6200 


0.0447 


150 


21.3395 


405.9141 


0.6121 


0.0492 


0.6103 


0.0483 


200 


18.2231 


355.5294 


0.6142 


0.0478 


0.6074 


0.0410 


250 


15.7638 


309.9113 


0.6122 


0.0481 


0.6075 


0.0388 


300 


14.7470 


295.3437 


0.5987 


0.0410 


0.5951 


0.0384 


350 


12.9422 


271.1457 


0.5865 


0.0350 


0.5845 


0.0386 


400 


12.1784 


249.9643 


0.5960 


0.0347 


0.5732 


0.0319 



TABLE I: Mean and standard deviation of probabihty of positive asymmetry with different 
fluctuation values, and mean of its single FRPs, as functions of N. Here, rj = 0.16, Ti = 200, Tr = 
20, and tq = 15. 

out of (dissipative) reversible dynamics. Furthermore, the weak dependence of our results on 
suggests that the asymmetries may survive in the large limit, with growing fluctuation 
values, although at present it is not possible to test directly very large A^'s and T's. If 
the asymmetries really survive, our results constitute a first (partial) verification of the 
predictions of Refs.^,0| (originally obtained for nonequilibrium stochastic systems), in the 



7 









P{a > 0) 


]E{a) 


250 


25.8539 


564.2185 


0.5840 


0.0305 


300 


23.8107 


558.0274 


0.6082 


0.0334 


350 


24.0250 


571.7431 


0.6219 


0.0484 


400 


20.3576 


554.6692 


0.6480 


0.0430 


350(*) 


40.6265 


1117.7722 


0.6297 


0.0474 


400(*) 


41.4094 


1099.1702 


0.6131 


0.0473 



TABLE II: Same as Table I, for ' only, with Te = 300. The cases with a (*) have Te = 500. 
context of time reversal invariant, deterministic dissipative particle systems. 
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